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(N . ABSTRACT 

We develop the mid-infrared extinction (MIREX) mapping technique of Butler & Tan (2009, Pa- 

■ per I), presenting a new method to correct for the Galactic foreground emission based on observed 
saturation in independent cores. Using Spitzer GLIMPSE 8 /im images, this allows us to accurately 
probe mass surface densities, S, up to :^ 0.5 g cm~^ with 2" resolution and mitigate one of the main 
sources of uncertainty associated with Galactic MIREX mapping. We then characterize the structure 
of 42 massive starless and early-stage cores and their surrounding clumps, selected from 10 infrared 
dark clouds (IRDCs), measuring ScK^) from the core/clump centers. We first assess the properties 
of the core/clump at a scale where the total enclosed mass as projected on the sky is Md = 60 Mq. 
We find these objects have a mean radius of Rd 0.1 pc, mean Sci = 0.3 g cm~^ and, if fit by a 

^ : power law density profile pd ex r ^^'^^ a mean value of kp^ci = 1.1. If we assume a core is embedded 

^ r"! ■ in each clump and subtract the surrounding clump envelope to derive the core properties, we find a 

O mean core density power law index of kp^c = 1-6. We repeat this analysis as a function of radius and 

I . derive the best-fitting power law plus uniform clump envelope model for each of the 42 core/clumps. 

2 ' The cores have typical masses of Mc ~ 100 Mq and Sc ^ 0.1 g cm~^, and are embedded in clumps 

with comparable mass surface densities. We also consider Bonnor-Ebert density models, but these 
^ . do not fit the observed S profiles as well as power laws. We conclude massive starless cores exist 

and are well-described by singular poly tropic spheres. Their relatively low values of H and the fact 
that they are IR dark may imply that their fragmentation is inhibited by magnetic fields rather than 
radiative heating. Comparing to massive star-forming cores and clumps, there is tentative evidence 
for an evolution towards higher densities and steeper density profiles as star formation proceeds. 
, Subject headings: ISM: clouds, dust, extinction — stars: formation 

m 

(N 

^ I 1. INTRODUCTION 

■ Characterizing the initial conditions of massive star formation is important for distinguishing between various the- 
' oretical models. For example, McKee & Tan (2002, 2003, hereafter MT03) have presented the "Turbulent Core 

Accretion Model", which assumes that massive stars (including binaries and other low-order multiple systems) form 
^ from massive starless cores in a manner that can be considered a scaled-up version of the standard theory of low-mass 
star formation (Shu, Adams, & Lizano 1987). The core undergoes global collapse to feed a central accretion disk. 
, These massive starless cores are also assumed to be near virial equilibrium and in approximate pressure equilibrium 
J-H ■ with the surrounding clump environment, whose pressure is likely to be set by the self-gravitating weight of the gas, 
, Pel GTi^Y-) where Ed is the mass surface density of the clump. The clump is defined to be the gas cloud that fragments 
to form a star cluster. Observed regions of massive star formation, including revealed massive star clusters where this 
activity occurred recently, have high values of Hci 1 g cm~^ (but with a dispersion of about a factor of 10), implying 
large values of Pd/k > 10^ K cm~^. 

A massive, virialized core in pressure equilibrium with this environment cannot be supported by thermal pressure, 
given observed temperatures of T ~ 10 — 20 K, and so must be supported by some form of nonthermal pressure, i.e. 
turbulence or magnetic fields. Since massive stars are "rare", in the sense that they constitute only a small mass 
fraction, ~ 5 — 10%, of the observed initial stellar mass function, then massive starless cores that will eventually form 
massive stars are also expected to be rare. Most massive structures are likely to fragment into clusters of lower- mass 
stars. We anticipate that preventing fragmentation of massive starless cores likely involves magnetic fields being strong 
enough such that the core mass is approximately equal to a magnetic critical mass and substructures are magnetically 
subcritical. Krumholz & McKee (2008) have argued fragmentation is prevented by radiative heating from surrounding 
low-mass protostars, which requires them to have high accretion rates and thus for the clump to have a high value of 
S > 1 g cm~^ in order to form massive stars. There is no such constraint if magnetic fields prevent fragmentation. 

MT03 modeled massive cores as singular polytropic spheres, with a power law density distribution of pc cx r~^P'''. 
There were few observational constraints on this density distribution, so MT03 assumed cores were part of a self-similar 
hierarchy of structure also shared by their surrounding clumps, where observations suggested kp^d — 1.5 (e.g. van der 
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Tak et al. 2000; Beuther et al. 2002; Mueller et al. 2002). The parameters of a fiducial massive core make it clear why 
it is difficult to measure the structure observationally. For a GOM© core embedded near the center of a clump with 

mean Sci = lgcm~^, the radius is Rc = 0.057M^^gQ(i;ci/l gcm~^)~^/^ pc. At typical distances, > 2 kpc, this radial 
size corresponds to < 5.9''. 

Until recently, measurements of mass surface densities of ~ 1 g cm~^ (equivalent to Nn = 4.27 x 10^^ cm~^ or 
Ay = 230 mag) were based mostly on mm dust continuum measurements, which require knowing the dust emissivity, 
dust temperature and dust-to-gas ratio. In particular, the dust emissivity and temperature may be expected to vary 
along the line of sight through a dense cloud. For total power observations with single dish telescopes, resolutions have 
been limited to ~ 10'', e.g. the diffraction limit of a 30 m telescope at 1.2 mm. Rathborne et al. (2006) carried out 
a study of this emission from 38 Galactic Infrared Dark Clouds (IRDCs) (these types of clouds are described in more 
detail below). The above considerations show that we do not expect their results to be able to resolve the scale of 
individual massive cores in high pressure environments, but are better suited to studying the properties of clumps that 
might form whole star clusters. Higher angular resolution has been achieved with interferometric observations, but 
these have been possible only towards relatively limited samples of objects, many of which are already in the process 
of forming a star (e.g. Bontemps et al. 2010). 

The advent of space-based, MIR, high-photometric-accuracy, imaging surveys of the Galactic plane has opened up 
a new way to probe high mass surface density structures. Indeed, these cold, high S structures were ffist identified 
as "Infrared Dark Clouds" from analysis of Infrared Space Observatory {ISO; Perault et al. 1996) and the Midcourse 
Space Experiment {MSX; Egan et al. 1998) imaging data. With the Spitzer Space Telescope^ more precise and higher 
angular resolution data became available. This led Butler & Tan (2009, hereafter BT09 or Paper I) to attempt to 
develop MIR extinction (MIREX) mapping as a precision technique for probing high mass surface density regions. 

BT09 adopted the Ossenkopf & Henning (1994) thin ice mantle coagulated (for 10^ yr at densities of 10^ cm~^ or 
equivalently for ~ 10^ yr at densities of ~ 10^ cm~^) dust model for their fiducial analysis and a gas-to- (refractory 
component)-dust mass ratio of 156 (slightly higher than the value of 141 estimated by Draine (2011) from depletion 
studies). With these values the opacity per unit gas mass in the IRAC band 4 at 8 /im of a source with a spectrum 
typical of Galactic diffuse MIR emission is hzsfim = 7.48 cm^ g~^. The fiducial value adopted by BT09 and in this paper 
is = 7.5 cm^ g~-^ so that 

TS/im 7.5 ( ) • (l) 

\gcm~^y 

From the variety of dust models considered by BT09, we expect ~ 30% uncertainties in the absolute value of hzsfim- 
Within a particular IRDC, we can expect some systematic variation in hcsfim due to different degrees of ice mantle 
growth, but these should be at most ~ 20% (Ossenkopf & Henning 1994), and probably much less after averaging over 
conditions on a line of sight through the cloud. 

Apart from the choice of MIR opacity, there are two main sources of uncertainty involved in MIR extinction map- 
ping. First, the intensity of the MIR emission behind the cloud is assumed to be smooth and must be estimated by 
extrapolation from nearby regions that are assumed to be extinction free. BT09 estimated that the extrapolation could 
lead to flux uncertainties of about 10%, corresponding to errors in S 0.01 g cm~^, given fiducial dust models. This 
is a minimum S, below which MIR extinction mapping becomes unreliable. These problems of background estimation 
can be reduced by choosing IRDCs that are in regions of the Galactic plane where the observed surrounding emission 
is relatively constant and smooth around the cloud. One systematic bias that we expect to be present is caused by the 
fact that there will typically be some cloud material in the "envelope" region around the IRDC where extinction was 
assumed to be zero. From studies of CO emission around IRDCs (Hernandez & Tan 2011; Hernandez et al. 2011) we 
estimate that this envelope typically has S 0.01 g cm~^. This is an additional reason why MIR extinction mapping 
becomes unreliable at low values of S. This problem can be addressed by combining MIR and NIR extinction mapping 
techniques (Kainulainen et al. 2011; Kainulainen & Tan, in prep.) The second major source of uncertainty is caused 
by foreground MIR emission along our line of sight to the IRDC. Neglecting this causes us to underestimate r and thus 
E. The effect can be minimized by choosing IRDCs that are relatively nearby, as was done by BT09. BT09 also tried 
to estimate the expected amount of foreground emission assuming it comes from a smoothly distributed population 
of small dust grains heated by massive stars that follow an exponential distribution in the Galaxy. For a cloud at a 
distance of 5 kpc at a Galactic longitude of / ^ 30°, we estimate that 27% of the observed Galactic diffuse emission is 
from material in front of the cloud. For a part of the cloud that has S estimated to be ~ 0.1 g cm~^ in the absence 
of a foreground correction, applying this correction raises S by about a factor of two. If this foreground correction is 
not applied then the largest values of E that can be derived are only ~ 0.2 g cm~^. This model-dependent estimate 
of the foreground is quite uncertain and one of the main reasons that BT09 concentrated on nearby IRDCs. Poor 
understanding of the foreground emission is likely to limit the reliability of the mass surface densities and masses of 
IRDCs derived for large samples of objects (e.g. Simon et al. 2006; Peretto & Fuller 2009), especially for the more 
distant objects. Local heating of small dust grains that then produce MIR foreground emission cannot be accounted 
for in the BT09 model of foreground estimation and this can affect even nearby IRDCs. However, choosing relatively 
quiescent IRDCs can help minimize this particular source of uncertainty. One of the main goals of this paper is to 
introduce a new, improved method to measure the intensity of the foreground emission. 

The structure of this paper is as follows. In ^ we introduce the method of saturation-based MIR extinction 
mapping. In ^ we present the results of applying this method to study the structure of 42 massive starless and 
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early-stage core/clumps located in 10 IRDCs. In §3]we discuss the implications of these results for massive star and 
star cluster formation theories. 

2. SATURATION-BASED MIR EXTINCTION MAPPING 

The MIREX mapping technique requires knowing the intensity of radiation directed towards the observer at a 
location just behind the cloud of interest, Ijy^o, and just in front of the cloud, /^^i. Then for negligible emission in the 
cloud and a simplified ID geometry, 

=e""^/.,o, (2) 

where the optical depth = hz^^, where is the total opacity at frequency z/ per unit gas mass and E is the gas 
mass surface density. 

We cannot see lu^o directly, so it must be estimated by interpolation from surrounding regions. BT09 tried two 
main ways to do this using median filters. The Large-Scale Median Filter (LMF) method used a square filter of size 
13' that was much larger than the IRDCs of interest so that the clouds did not significantly depress the estimated 
median intensity. This has the advantage of not assuming any prior knowledge about the IRDC, but the disadvantage 
of a coarse angular resolution of background intensity fiuctuations. For studying specific IRDCs that can be defined 
as occupying a certain region of the sky, e.g. an ellipse, we thus introduced the Small-Scale Median Filter (SMF) 
method. Here the size of the filter is set to be one-third of the major axis of the IRDC ellipse (defined by Simon et 
al. 2006), but it is only applied for background estimation outside of the IRDC ellipse. Inside the ellipse we estimate 
the background by interpolating from the surrounding background model. BT09 estimated that the uncertainties in 
background estimation due to this interpolation were at a level of < 10%, which corresponds to S < 0.013 g cm~^. 

However, because of foreground emission towards the IRDC, we actually observe (see Figure [T]) 

-^z^,l,obs — Iu,fore ~^ — -^i/,fore ~^ 6 7 (3) 

and towards the IRDC surroundings, where we are trying to estimate /^y^o, we actually observe 

^zy,0,obs = Iu,fore + lu,0' (4) 

The primary uncertainty in the MIREX mapping method of BT09 for larger values of S is the estimate of the level of 
the foreground contribution to the intensity, Iu,iore- In order to increase the method's sensitivity to higher values of 
mass surface density, we now describe a new, empirical method to estimate this contribution. 

If a core has a high enough mass surface density, then it will block essentially all the background emission. The 
observed minimum intensity in the cloud will then be approximately equal to the foreground emission and the angular 
distribution of this intensity may appear to fiatten or "saturate" . It is difficult to be certain if an individual dense core 
is saturated (as is sometimes assumed if the foreground is simply estimated from the darkest pixel, e.g. Ragan et al. 
2009). However, we propose that if the minimum intensity is observed to "be the same" in two or more "independent" 
cores (i.e. spatially resolved peaks in S) in the same cloud, then this is very likely due to saturation. In practice, by 
"be the same" we adopt the condition to be within 2a of each other, where a is the uncertainty in the GLIMPSE 8 /im 
intensities of 0.6 MJy sr~^ (Reach et al. 2006). By "independent" we adopt an angular separation of at least 8'' i.e. 
much larger than the 2'' FWHM of the Spitzer IRAC 8 /im PSF. 

The algorithm for this method is as follows (see also Figure [T]): 

1) Define a region of the sky as the "IRDC". Following BT09, we use the ellipses from the catalog of Simon et al. 
(2006), which were based on MSX images. 

2) Using GLIMPSE 8 /im images, find the minimum value of Iu,i,ohs inside the IRDC, /iy,i,obs(min). 

3) Search for all pixels in the IRDC with /iy^i^obs(niin) < /zy,i,obs < ^z/,i,obs(niin) + 2a. If there are pixels meeting this 
criteria that are independent (to be conservative we adopt > 8'' away from the IRDC minimum), then the IRDC is 
defined to be saturated, all the above pixels are labeled as "saturated pixels", and the following steps are carried out. 

4) The mean value of /z/,i,obs of the saturated pixels is evaluated, /i/,i,obs(sat). We set the foreground intensity 
(which includes all sources of emission: Galactic, Zodiacal and instrumental, Battersby et al. 2010) to be Iu,fore = 
^zy,i,obs(sat) — 2a. This subtraction is motivated to have Iu,fore < ^zy,i,obs(min) and thus give every pixel a finite value 
of r and thus S. 

If all the "saturated pixels", defined above, really did have negligible values of Ij^^i and had a distribution of 
intensities that was relatively uniform in the above range, then Iv^fove — ^z/,o,obs(niin) + la and our method would be 
underestimating Iv.ioie by 2a = 1.2 MJy/sr. In fact, we do find for the ^ 300 "saturated pixels" in the 10 IRDCs of 
our sample, a mean value of /zy,i,obs — ^zy,i,obs(min) 0.7 MJy/sr. Thus, we are likely to be underestimating lu^fore 
(overestimating /j,^i) by an amount :^ 2a = 1.2 MJy/sr. This leads to a value of S where our measured values are 
significantly affected by saturation: 

where h,i/Iv,Q = Vm,i/Vm,o = e-^^^-^"^*) so r8^m(sat) = ln(/^,o/^z.,i), that is set by = 2a ^ 1.2 MJy/sr. 
For a typical IRDC with /r/,o,obs = 100 MJy/sr, Iu,iore = /fore^i/,o,obs = 30 MJy/sr so that /j,,o = 70 MJy/sr, then 
T8/im(sat) = 4.07 and S(sat) = 0.544 g cm~^. For a region of such a cloud with a true value of S = 0.5 g cm~^ so that 
/jy^i = 1.65 MJy/sr, if we have underestimated /zy,fore by 1.2 MJy/sr, then we would infer S = 0.427gcm~^. Similarly, 
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for a true S = 0.4gcm~^, we would infer S = 0.361 gcm~^. The values of i;(sat) calculated with /g/im,! = 1-2 MJy/sr 
for the 10 IRDCs in our sample are listed in Table 1. They range from S (sat) =0.33 to 0.52 g cm~^ as one progresses 
along the Galactic plane towards / = 0, where the background is brightest. 

An additional uncertainty results from our use of a single effective value of /^g/xm = 7.5 cm^ g~^ averaged over the 
Spitzer IRAC 8/im band, weighting by the filter response function, the spectrum of the Galactic background and the 
dust opacity model (BT09). Since these functions vary over this wavelength range (see Fig. 1 of BT09), at large optical 
depths the actual transmitted intensity will be greater than that predicted, being more dominated by the region of 
the spectrum with the lowest opacity. The net effect is an underestimation of the true mass surface density, given 
the observed ratio of transmitted to incident intensities. We have investigated the size of this effect by integrating 
the transfer equation ([2]) over the above weighting functions (see Fig. [2]). For our fiducial dust model (the moderately 
coagulated thin ice mantle model of OH94), which has a relatively flat MIR opacity law, the effect is small: just a 
few percent effect up to value of S ~ 1 g cm~^, rising to about a 10% effect by S = 10 g cm~^. For illustrative 
purposes. Fig. [2] also shows the results for the Draine (2003) Rv =3.1 dust model, more appropriate for the diffuse 
ISM, which has bare grains and stronger variation of opacity across this wavelength range. Now the effect leads to an 
underestimation of S by up to several tens of percent for H ~ 1 g cm~^. Other dust models we have considered, such 
as the Draine (2003) Ry = 5.5 model, have somewhat smaller underestimation factors. 

A further additional systematic uncertainty results from the fact that the foreground intensity will vary across the 
IRDC, especially due to local radiation sources. The accuracy of the S values will be higher in regions closer to the 
locations of saturated cores, where Iu,iore has been estimated and in IRDCs with minimal local heating sources. 

3. RESULTS 
3.1. IRDC Properties 

Following the above algorithm, we find that all 10 IRDCs of the BT09 sample exhibit the effects of saturation. 
In hindsight, this is not too surprising since these clouds were selected to have relatively high contrast against the 
background. The S maps of the clouds are shown in Figs. [3] and |H The properties of these clouds are listed in Table 
1, where we also compare their properties to those derived with the SMF method of BT09 with the analytic model of 
foreground estimation. Using the saturation-based estimate of foreground emission, we find Iu,iore and thus /fore has 
increased in all the clouds. Thus the highest values of S that we infer have risen from ~ 0.1 — 0.3 g cm~^ in BT09 to 
~ 0.4 — 0.6 g cm~^ in this paper. The mean values, Ssmf rise by smaller factors, so that the total cloud masses rise 
by on average a factor of 2.0. A comparison of the global properties of these IRDC with the predictions of theoretical 
models of the interstellar medium will be presented in a separate paper. 

3.2. Massive Starless Cores and Clumps 
3.2.1. Locating the Cores 

The cores we are considering are a subset of those originally identified by Rathborne et al. (2006) based on their mm 
dust continuum emission, observed with the IRAM-30m Telescope at 11'' FWHM angular resolution. BT09 selected 
43 cores from the Rathborne et al. sample, excluding those with significant 8 /im emission and those with low-contrast 
against the MIR background (i.e. with E < 0.02 gcm~^). Here, we have excluded one of the BT09 cores, E4, because 
its GLIMPSE image suffers from a diagonal boundary artifact where the intensity of the diffuse emission changes 
abruptly. 

BT09 treated the cores as circular with radii equal to half the reported FWHM diameter of Gaussian fits that 
Rathborne et al. (2006) fitted to their mm continuum images. These circles were centered at the coordinates estimated 
by Rathborne et al. As discussed in ^3.1\ we expect our derived values of S to be higher (and more accurate) than 
those of BT09. Comparing the core masses of BT09 with those derived here for the same regions, we find they have 
typically increased by a factor of about 2.2. 

In this paper, we now redefine the core center to be the center of the highest S pixel inside the previous core 
boundary. If there are two or more adjacent saturated pixels at the core center, then their average position is used to 
define the center. In fact, 17 of the 42 cores exhibit saturation. Occasionally, after inspecting the 8 jam GLIMPSE and 
24 /im MIPSGAL images, we note the presence of MIR sources near (< 7.5") the core center. This occurs in 9 of the 
42 cores (B2, C6, C8, D5, D6, D8, E2, E3, II). In order to focus on massive starless and early-stage cores, we shift 
the center to a new, nearby (< 3'') E maximum to avoid any major sources of MIR emission within a radius of 7.5'' 
of the new center. In several cases (C4, D4, F2, Jl), the E map inside the Rathborne et al. core boundary does not 
exhibit a well-defined high S peak. In these cases we select a new core center as close as possible to the Rathborne et 
al. core: normally this is within a few arcseconds of the boundary, but for F2 it is about 10" outside. 

Figure [5^ shows the S map of core Al, extracted from the larger image of IRDC A, shown in Fig.[3K. Pixels suffering 
from saturation are marked with small white squares. The core center is marked with a cross. Similar images of all 
42 cores are shown in Figs. [6lto[T2l 

We note that 5 of the IRDCs (B, E, G, H, J) only have one core that exhibits saturation. This is possible because 
the condition to determine if an IRDC is saturated is based on independent positions (separated by at least 8") having 
the same foreground intensity (to within 2cr), rather than requiring 2 cores to meet this condition. The 42 cores we 
have selected for analysis are not meant to be a complete census of all the dense regions in these IRDCs. For example, 
IRDC J only has one core selected. 
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The core S maps exhibit complex structure. It is not easy to define the boundary of a core from its surrounding 
clump, especially when one recalls we are viewing a 3D structure in projection. CO (1-0) data exist for these IRDCs 
via the Galactic Ring Survey (Jackson et al. 2006), but with poor angular resolution 46''). Also, in the cores 
we expect CO to be highly depleted from the gas phase due to freeze-out onto dust grain surfaces. Widespread CO 
depletion has been observed in IRDC H by comparing our S map with C^^O emission observed with the IRAM-30m 
Telescope (Hernandez et al. 2011). Fontani et al. (2011) observed N2H+, which does not freeze-out so readily as CO, 
from 4 of our cores (Cl, Fl, F2, G2), but again with relatively poor angular resolution (> 10''). Thus, given the lack 
of high angular resolution molecular line data for all the cores, here we present a uniform analysis of core structure, 
based only on the extinction maps. 

For simplicity, we first make radial profiles of mean total mass surface density, which we refer to as Sd since it 
includes contribution from the clump (see below), considering a series of annuli extending from the core center with 
width equal to 1 pixel, i.e. 1.2" . Fractional overlap of pixels with these annuli are accounted for. "Holes" in the 
S maps due to MIR sources are treated as having a zero, i.e. negligible, value. In general, these sources do not 
significantly affect our characterization of core structure, at least in the inner ~ 7.5", since we have chosen cores that 
are relatively free of strong sources (E3 is the worst affected, and is somewhat exceptional in this regard). Larger 
annuli are minimally affected by individual MIR sources, which cover only a small fraction of the area. We extend the 
radial profiles out to a maximum angular scale equal to that reported by Rathborne et al. (2006) based on mm dust 
emission, i.e. a radius equal to one FWHM of their fitted Gaussian profile. As we will see, this is generally larger than 
the scale over which the core can be considered to be a single monolithic object. 

For Core Al, Figure [SJd shows Sci(^) with blue open square symbols, plotted at the radii corresponding to the 
center of each annulus. The total enclosed mass, which we refer to as the clump mass Mci(r), is indicated by the blue 
long-dashed line. 



Before considering a more detailed analysis of the radial structure, it is instructive to first consider the properties 
of these core/clump objects at a scale where the total mass enclosed is M^i = 60 Mq. If all this mass were in a core, 
then such a core has the potential to form a ~ 30 Mq star, given expected star formation efficiencies of ~ 50% due 
to protostellar outflows (Tan & McKee, in prep.). Note that there is of course no guarantee that all our sources will 
collapse in this way and on statistical grounds one would not expect them to: most are likely to undergo fragmentation 
to form lower-mass starsQ For Core Al, 60 Mq is enclosed within Rd = 0.0962 pc, so at this scale Ed = 0.431 gcm~^. 
A black dashed circle with this radius is shown in Fig. [5^ and these core properties are listed in Table 2. Core Al 
happens to be one of the most extensively saturated cores at the scale of an enclosed mass of 60 Mq (along with C2, 
HI, II, 12, Jl), so these numbers are likely to be significantly affected by saturation (which causes us to underestimate 
S), so actually the radius enclosing 60 Mq would be smaller and S larger. 

The distributions of the radii, Rd, and mean mass surface densities, Sci, of the 42 core/clumps at the M^ = 60 Mq 
scale are shown in Figure [T3k (with the 6 highly saturated cores — Al, C2, HI, II, 12, Jl — shown as a shaded 
subset). The mean/median/RMS dispersion-about-the-mean of Rci{Md = 60 Mq) = 0.121/0.114/0.0238 pc. The 
mean/median/RMS dispersion values of Sci(Mci = 60 Mq) = 0.296/0.318/0.0952 g cm-^ (see also Table 3). 

We next fit a power law density distribution. 



where ps,ci = Mh^h,s,c1 (with //h = 2.34 x 10~^^ g) is the density at the surface of the clump, Rd. We project the 
above distribution to derive Sci(r), which we then convolve with a Gaussian with a FWHM of 2" (to allow for the 
Spitzer IRAC 8 /im PSF). We then fit this model to the observed Sci(^) profile, excluding annuli that are significantly 
(> 50%) affected by saturated pixels. For Core Al, kp^d = 1-40 and nH,s,ci = 2.47 x 10^ cm~^. For the whole sample, 
the mean/median/dispersion values of kp^d = 1.09/1.10/0.236 and nH,s,ci = (1.76/1.85/0.852) x 10^ cm~^. These 
distributions are shown in Fig. [T3k with the blue dotted histograms. The values for individual cores are listed in Table 
2. 

The above analysis is somewhat simplistic in that it has assumed the structure exists in isolation. In reality, we see 
that these high S objects are surrounded by regions that also have significant mass surface densities. Thus next we 
model the cores with a similar power law density structure. 



but now when comparing to the observed S maps we account for the mass surface density of the surrounding clump 
medium, Slci,env We estimate Sci,env using the observed value in the annular region from Rc to 2Rc. This choice is 
motivated by the desire to sample a region of the clump that has a scale comparable to the core in both size and 

^ One cannot distinguish between massive star formation Core Accretion and Competitive Accretion theories (Bonnell et al. 2001) simply 
by observing that a massive structure is actually composed of sub- fragments (c.f. Bontemps et al. 2010, their section 4.6). 



3.2.2. Core and Clump Properties at the 60 Mq Enclosed Mass Scale 




(6) 




(7) 
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masfl We assume this same value of i;ci,env covers the area of the core, and so subtract it from the interior T^cii^) 
profile to derive the mass surface density profile of the core, Hc(^)- Thus note that Sci,env = — ^c- 

For the same Md = 60 Mq enclosed mass scale as defined above, we set Rc = Rci{Mci = 60 Mq). The masses 
contained in the cores, based on integrating the Hc(^) profile, are of course less than the 60 we previously identified 
with the clump. For Core Al, we derive a core mass of Mc = 37.9 with Sc = 0.316gcm~^. Its value of kp^c = 2.04 
and nH,s,c = 1.21 x 10^ cm~^, i.e. a steeper density profile with a lower value of the volume density at the surface 
than was derived previously. For the 42 cores we find the mean/median/dispersion values of Mc = 30.4/30.4/12.7 M©, 
Sc = 0.139/0.160/0.0738 gcm-2, kp^c = 1.64/1.67/0.271 and nH,s,c = (0.639/0.750/0.394) x 10^ cm-^ (See red solid 
line histograms in Fig. [T3k and Tables 2 & 3). Compared to the clump results (i.e. derived from the total S profiles), 
above, for the envelope-subtracted core properties we necessarily find smaller surface densities, steeper density profiles 
and smaller volume densities. 

The power law fits ignore annuli affected by significant saturation, where S is underestimated. Thus we also estimate 
a core mass, Mc,pl, based on extrapolation of the power law fits to the center of the core: 

There are no cores where the derived kp > 3 for which the inner boundary condition would have be considered. If 
there were, then in these cases we would expect to truncate the power law at the Jeans scale in the core. For Core Al, 
^c,PL = 49.2M0, about 30% times larger than Mc. Such an increase is expected since this is one of the most extensively 
saturated cores. For the rest of the 42 cores the change is typically much smaller. The mean/median/dispersion values 
of Mc PL = 31.1/31.0/13.5 Mq and Sc pl = 0.154/0.171/0.0899 (See red solid line histograms in Fig. ^ and Tables 
2 & 3). 

3.2.3. Best-fit Power Law Cores 

We now repeat the power law core plus clump envelope fitting procedure as a function of radius, starting at the 
inner region with 3 unsaturated annuli for the core. An annulus twice as large in radius is used to estimate the value 
of S of the clump envelope, Sci,env We assess the relative goodness of fit of this model as a function of r by finding 
the minimum of the reduced parameter, defined by 

V- 1 Fc,PL(r) - ^cAr)? .Q^ 

where N is the number of annuli, u = N — 2 is the number of degrees of freedom and a is the error for each annulus, 
which we take to be a = 0.01gcm~^ + 0.2i;c. Note, that because of the 2'^ angular resolution of Spitzer IRAC^ adjacent 
annuh are not completely independent. However, the relative values of should still give a measure of the best-fitting 
model. 

We place some additional constraints on the fitting. First, we do not allow the best-fit core to extend beyond 
neighboring core centers (from our sample of 42 cores). Second, to prevent independent discrete structures that are 
not part of our core sample from influencing the fitting, we check for a 3<j rise in the profile by comparing the 
annulus before any rise begins to the following local maximum. If this occurs, we ignore fits beyond the pre-rise 
annulus, and search inward for a local maximum in xi ^ind define that to be the best-fit radius (for example, this 
occurs in Core A3). Third, if more than 25% of an annulus is composed of MIR emission pixels, we do not extend the 
fit any further. In these cases, a prior unaffected annulus with a local maximum in xt is chosen as the best-fit radius 
(this circumstance only arises in Core E3). As a result of the above constraints, it is possible that the global minimum 
of xt will not be chosen as the "best-fit" . 

The results for Mc(r), kp^d^) and — logx^(^) are shown for Core Al in Fig. [SJd. The location of the peak value 
of — indicates the best-fitting power law (PL) core radius, which occurs at 0.251 pc with a value of = 1.62. A 
circle of this best-fit core radius is shown in the S map of the core in Fig. [5^. The total enclosed mass at this scale 
is Mci = 303 M0, the core mass is Mc — 194 M^, the mean core mass surface density is = 0.204 g cm~^ and the 
clump surrounding the core has i;ci,env = ^ci — = 0.115 gcm~^. The core mass based on integrating the power law 
profile is Mc,pl = 204 M©, yielding a slightly higher mean mass surface density of Sc,pl = 0.214 g cm~^. 

The best-fit total Sci(^) = ^c(^) + Sci,env model profile is shown by the solid line in Fig.[5)3 (the dotted continuation 
in the inner region indicates where annuli affected by saturation are not used in the fitting) . Figure [5b shows the 
clump envelope subtracted profile of Sc(r), together with various projected power law fits, including the best-fit value 
of kp^c = 1-88. The parameters of the best-fitting power law plus clump envelope model are listed in Table 2. 

The distributions of i?c, Sd (which is the mean total S over the area of the core), Sc,pl, ^p,ci5 ^p,c7 ^H,s,cb ^h,s,c5 
Mci and Mc,pl are shown in Figure [T3b and summarized in Table 3. The values for each core are listed in Table 2. 

It is important to note that these "best-fit" values may not necessarily be the most accurate description of the core 
structures. They are based on azimuthally- aver aged quantities. The H map of a particular core should be inspected 

^ We have also tried measuring Ed env from a thin, 1.2^^ wide annulus just outside which generally leads to larger estimated values 
of Eh 

env However, we consider that this thin-shell annulus does not sample a large enough region and mass of the clump that, via 
self-gravity, would be responsible for setting core's surrounding pressure (see discussion in ^T\i. 
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to gauge the validity of this assumption. Also, the values of as a function of radius should be checked to gauge the 
reasonableness and uniqueness of the fit. 

The radii and masses of the best-fit cores are generally, but not always, larger than those at the M^i = 60 Mq scale, 
and thus the volume densities are generally lower. The mean/median/dispersion values of kp^c\ = 1.10/1.12/0.246 and 
kp^c = 1.58/1.56/0.277 are however very similar to those derived at the M^i = 60 Mq scale, which suggests that the 
assumption by McKee & Tan (2002) and MT03 of a self-similar hierarchy of structure from the clumps to core scales is 
a reasonable one. The fiducial value they adopted of kp = 1.5 also is close to the average values found in this sample. 
In Figure [HI we plot kp^c versus Mc,pl, ^c,pl and Hci,env for the best-fit cores. There are no apparent correlations of 
kp^c with these properties. 

3.2.4. Best-fit Bonnor-Ebert Cores 

We perform a similar analysis as the power law plus constant envelope fitting as a function of radius, but now using 
critical Bonnor-Ebert profiles (varying the total effective sound speed, and surface pressure, Pq) plus a constant 
envelope. See Dapp & Basu (2009) for more details about fitting Bonnor-Ebert profiles to column density data. 

The best-fitting model for Core Al is shown in Fig.[5li. This has Rc = 0.670 pc, = 353 M©, Se = 0.0523 gcm-^ 
Cs = 0.275 km s~^ and Po/k = 8.9 x 10^ Kcm~^. However, the value of = 8.75, which is significantly larger, 
i.e. worse, than the best-fit power law plus clump envelope model fit (for which = 1.62). Also the size of the 
Bonnor-Ebert fitted core is much larger than the power law model: as can be see from Fig. [5^, on these larger scales 
the assumption of single monolithic and azimuthally symmetric structure becomes less valid. 

Carrying out the Bonnor-Ebert analysis for all 42 cores, we find the fits are generally worse than for the power law 
models. The best-fitting Bonnor-Ebert radii are typically larger than those of the power law core models. For these 
reasons, we do not consider the Bonnor-Ebert models further in our discussion. 

4. DISCUSSION & CONCLUSIONS 

We have presented a new method to accurately probe mass surface densities in the range ~ 0.01 to ~ 0.5gcm~^ on 
arcsecond scales in quiescent, infrared dark clouds, some of which are likely to be the sites of future star formation. The 
method uses the small-scale median filter method of background interpolation from regions around a defined IRDC 
(BT09) and then estimates the level of foreground emission by seeing if there are independent, nearby, saturated cores 
within the IRDC. If so, the foreground level is set equal to that observed towards these saturated regions. The resulting 

5 measurements derived from this MIREX mapping depend on the assumed MIR dust opacity per unit total mass, 
but do not depend on the dust temperature, which is a distinct advantage over measurements based on sub-mm/mm 
dust continuum emission. 

Focusing on 42 core/clumps within 10 IRDCs, we have tried various methods of characterizing their azimuthally- 
averaged structure. Our preferred method, following the model of McKee & Tan (2002, 2003), involves fitting power 
law cores surrounded by a clump envelope, which is assumed to have a constant value of S that can be estimated from 
the surrounding region. We have fitted these models as a function of radius from the core center, deriving an overall 
best-fit, but also presenting the full results of this radial characterization. The typical value of the volume density 
power law index that best describes the cores is kp^c — 1-6. This is close to the fiducal value of 1.5 adopted by McKee 

6 Tan (2002, 2003), who based their choice on previous measurements on the larger, ^parsec, scales of gas clumps. 
We find this power law index does not appear to vary significantly with scale (i.e. between the 60 Mq enclosed mass 
scale and the best-fit power law plus clump scale), nor with other core or clump properties, suggesting the presence of 
a self-similar hierarchy of structure. 

On the scale at which the total projected enclosed mass is 60 M©, the derived cores have about 50% of this mass. 
If massive star formation is to occur, then this is the material that has a high probability of being incorporated into 
the massive star. These cores have typical radii of 0.1 pc, masses of ~ 30 M©, mean mass surface densities of 
He ^ 0.15 gcm~^ and surrounding clump mass surface densities of similar values. If one regards our method of clump 
envelope subtraction to be an overestimate, then one can consider the typical properties of the clumps on these scales 
as being representative of the gas that will form massive stars, i.e., with Md = 60 Mq and — 0.3 g cm~^. 

The above values of Hci are lower by factors of ~ 3 — 7 than the fiducial value of 1 g cm~^ considered by McKee 
& Tan (2002, 2003). Note, their theoretical model is general and does not require a particular value of Ed, so this 
difference does not require any physical explanation. However, Krumholz & McKee (2008) have proposed massive star 
formation requires ^ lgcm~^, based on a model in which fragmentation of massive cores is prevented by radiative 
heating from surrounding lower-mass protostars. The high value of Sd is required so that the lower-mass protostars 
accrete at high enough rates that they are luminous enough to sufficiently heat the massive core. 

Since massive star formation occurs relatively rarely, it may be that it occurs preferentially in cores with higher values 
of Sc and Hd than we have observed for the average of our sample, which we note does show significant dispersion. 
However, we consider this unlikely given that we have selected the highest S regions from 10 IRDCs that show some 
of the highest contrast against the Galactic MIR background (selected from the larger sample of 38 IRDCs studied by 
Rathborne et al. 2006). 

Saturation limits our S maps to values of c± 0.5 g cm~^ and this may be affecting our ability to find the highest S 
cores. However, we are excluding the saturated regions when deriving core and clump density profiles and most cores 
do not exhibit extensively saturated centers. We expect our choice of MIR opacity per unit gas mass may be uncertain 
by ^ 30% (see §1), so this by itself is unlikely to explain the relatively low values of S that we are deriving compared 
to the Krumholz & McKee (2008) prediction. 
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Another possibility, is that these cores and clumps will evolve to higher values of S before massive star formation 
occurs. Indeed, the fact that these are IR dark objects suggests that they cannot yet be experiencing much radiative 
heating. However, there is observational evidence for star formation activity in some of these cores. For example, Y. 
Wang et al. (2006) observed water maser emission located 4.3f from the center of Core CI. K. Wang et al. (2011) 
have reported protostellar outflows from Core C2. 

To better compare our IRDC core/clump sample with more evolved systems, in Figure [T5l we show the 42 core/clumps 
on the S versus M diagram (following Tan 2007). Here, S is measured from the total observed mass inside a given 
radial distance from the core/clump center. In Figure [TBI we compare these profiles to the properties of the 31 star- 
forming clumps whose IR and sub- mm dust continuum emission was observed and modeled by Mueller et al. (2002). 
Note that these properties depend on the (ID) modeled temperature structure, dust emissivity (they used the same 
Ossenkopf & Henning (1994) dust model that we have adopted for our MIREX maps) and gas-to-dust ratio (we have 
scaled Mueller et al.'s masses by a factor 1.56 to be consistent with our adopted gas-to-dust ratio). 

The IRDC cores/clumps overlap only with the lower-S range of the star- forming core/clump sample, perhaps indi- 
cating there is a (physically plausible) evolutionary growth in core/clump density as star formation proceeds. However, 
note that many (indeed most) star- forming cores and clumps have S < 1 gcm~^. Alternatively, the lack of starless 
high H core/clumps may be due to the somewhat smaller volume of the Galaxy that we have probed with our nearby 
IRDC sample, compared to the Mueller et al. star- forming core/clump sample. 

Mueller et al. (2002) found density power law indices of kp^c\ = 1.8 ± 0.4, slightly steeper than our derived values for 
IRDC cores of kp^c — 1-6, but significantly steeper than our value for clumps of kp^^i — 1.1. Again, this latter difference 
may indicate an evolution in cloud properties as star formation proceeds. 

Considering the above results, we suggest that the initial conditions of local massive star formation in the Galaxy 
may be better characterized with values of Sd — 0.2gcm~^ rather than lgcm~^, which would imply smaller accretion 
rates and longer formation times that the fiducial values of MT03. The accretion rate becomes 

= 1.37 X 10-^ f ^ V'V Meyr- (10) 

for a core with kp^c = 1.5 and a star formation efficiency of 50%, where is the instantaneous protostellar mass and 
m*/ is the final protostellar mass. The star formation timescale becomes 

t*.=4.31xl0M^^) ( A yr. (11) 

In this case of massive star formation at relatively low values of S, we expect fragmentation of the cores is prevented 
by magnetic fields, i.e. if the core mass is equal to the magnetic critical mass (Bertoldi & McKee 1992) 

where R and Z are the major and minor axes of the core, B is the mean field strength in the core, and nn is the 
mean number density of H nuclei. Thus for a core with nn = 10^ cm~^ and Mc = 100 M©, typical of our sample, the 
condition Mb = Mc requires a field strength 



100 Mo y \Rj VlO^ 

If cores have some significant magnetic support, then we expect R/Z > 1, perhaps ~ 2, so that the required field 
strength in eq. (p!3|) is then 190 jaG. Such field strengths are similar to those observed in regions of active massive star 
formation (e.g. Crutcher 2005). Indeed, Crutcher (2005) noted the observed mass to flux ratios scattered about the 
critical value. Numerical simulations of the collapse of marginally magnetically critical (rather than super critical, e.g. 
Wang et al. 2010; Hennebelle et al. 2011) cores are required to investigate this scenario for forming massive stars, 
and, more generally, for explaining the high-mass tail of the initial mass function (Kunz & Mouschovias 2009). 
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Fig. 1. — Schematic of simple ID model of radiative transfer through an IRDC, assuming negligible emission from the IRDC at frequency 
u. If independent cores (i.e. localized density maxima) A and B are both of sufficiently high E, then I^^^i <C lu^iore — ^i/,i,obs(^5 
providing an accurate, empirical estimate of the foreground intensity to the IRDC. 
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Fig. 2. — Effect of finite filter width on estimates of E and accuracy of an approximation using a single band-averaged opacity (see text). 
The Spitzer IRAC band 4, i.e. 8/im, filter has sensitivity from about 6.5 to 9.5 /im. Top panel: Ratio of transmitted to incident flux as a 
function of true mass surface density, Strue- The result for the band- average opacity for the moderately coagulated thin ice mantle dust 
model of OH94 (our fiducial model) is shown by the lower red dot-dashed line. The actual transmitted flux, calculated by integrating the 
transfer equation over the bandpass, is shown by the upper red dot-dashed line. The equivalent quantities for the Draine (2003) Ry =3.1 
dust model are shown by the lower and upper black solid lines: the effect is larger here as this dust model shows larger opacity variations 
across the band. Bottom panel: Effect on estimation of E. Given an observed ratio of transmitted to incident intensities, the true mass 
surface density, Strue, will be greater than that estimated using the band- averaged opacity, E. The ratio of Etrue/^ is shown by the red 
dot-dashed line for the above OH94 thin ice mantle model. The OH94 uncoagulated thin ice mantle model gives essentially the same result. 
The error is a few percent in the region of interest of the IRDC cores in this study. Also shown are these effects for the Draine (2003) 
Rv =3.1 (black solid line), Ry = 5.5 (black dotted line) and Ry = 5.5 Case B (black dashed line). 
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Fig. 3. — Mass surface density, Ssmf, maps of IRDCs A-F derived from MIREX mapping using Spitzer IRAC 8 /xm images with pixel 
scale of 1.2^^ and angular resolution of 2^^ using a saturation-based estimate of the foreground emission (^. The color scale is indicated 
in gcm~^. The dashed ellipse, defined by Simon et al. (2006) based on MSX images, defines the region where the background emission 
is estimated not directly from the small-scale median filter average of the image intensity, but rather by interpolation from nearby regions 
just outside the ellipse. The locations of the massive starless cores we have selected for analysis are marked with crosses. Bright MIR 
sources appear as artificial "holes" in the map, where we have set the values of E = g cm~^. 
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Fig. 5. — (a) Top Left: Mass surface density, Ssmf, map in gcm~^ of IRDC Core Al, extracted from the map of IRDC A (Fig. O. The 
core center is marked with a cross. Saturated pixels, for which E is a lower-hmit of the true value, are marked with small white squares. 
The black dashed circle shows the radius enclosing a total mass of 60 Mq . The red solid circle shows the extent of the core derived from the 
best-fit power law (PL) core plus envelope model (see text), (h) Bottom Left: Radial profiles of Core Al: observed log Eci/(gcm~2) (blue 
open squares, plotted at annuli centers) derived from the map shown in (a); total projected enclosed mass, M^i, (blue long-dashed line [see 
right axis]); core mass, Mc after clump envelope subtraction (red dashed line [see right axis]); index of core PL density profile, kp^c-, (red 
crosses); — logx^ (red triangles) of the PL plus envelope fit (best-t has a maximum or local maximum value [see text]); the best-fit PL plus 
envelope model (blue solid line; dotted line shows range affected by saturation that was not used in the fitting); log Ec/(gcm~^) of best-fit 
core after envelope subtraction (red solid squares) and PL fit (red solid line; dotted line shows range affected by saturation that was not 
used in the fitting), (c) Top Right: Sc(r), i.e. after clump envelope subtraction for the best-fit model (red solid squares; open squares show 
residual, post-subtraction envelope material). PL models with various values of kp^c are indicated (dashed lines), including the best-fit 
model with kp = 1.88 (solid line), (d) Bottom Right: As for (c) but for Bonnor-Ebert (BE) plus envelope fitting. ^c(r), i.e. after clump 
envelope subtraction for the best-fit model (red solid squares). Best-fit BE model (solid line) and models varying Cg (long-dashed lines) 
and Pq (dashed lines) by factors of 2 from this are shown (see text). 




Fig. 6. — Core Al, A2, A3, Bl, B2, CI E maps (notation as Fig. [St-) and azimuthally averaged radial profile figures (notation as Fig.fSj^). 




Fig. 7. — Core C2, C3, C4, C5, C6, C7 E maps (notation as Fig.[5ji) and azimuthally averaged radial profile figures (notation as Fig. [St). 
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Fig. 8. — Core C8, C9, Dl, D2, D3, D4 E maps (notation as Fig.[5ji) and azimuthally averaged radial profile figures (notation as Fig.[5j)). 




Fig. 9. — Core D5, D6, D7, D8, D9, El E maps (notation as Fig.[5ji) and azimuthally averaged radial profile figures (notation as Fig.lSj)). 
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Fig. 13. — (a) Left column: Distributions of properties of clumps and cores at the scale where the observed total enclosed mass Md = 
60 M0, which defines Rd = Rc- The first graph shows the distribution of Rd (red solid line) of the 42 cores enclosing this mass. The shaded 
subset shows the 5 cores (Al, C2, HI, II, 12) that have extended saturation in their centers on this scale. The second panel down shows 
the mean mass surface density of the clumps, Ed (blue dotted line; shaded subset as before), and cores after clump envelope subtraction, 
^c,PL (red solid line), based on the power law fit. The third panel down shows the distributions of kp^cl (blue dotted line) and kp^c (red 
solid line). The fourth panel shows the distributions of nH,s,cl (blue dotted line) and nH,s,c (red solid line). The bottom panel shows the 
distribution of Mc,pl (red solid line). The vertical blue dotted line shows the 60 Mq scale of the clump, (b) Right column: As for (a) but 
now for the best-fit power law plus clump envelope models. The only difference is that in the bottom panel, the clump mass, M^i, is now 
shown (blue dotted line). 
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Fig. 14. — /cp,c versus Mc^pl, Sc,pl and Sci,env for the best-fit cores. 
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Fig. 15. — Mass surface density versus mass (E — M) diagram (including lines of constant radial size and density [nu] for spherical 
clouds) for the 42 IRDC core/clumps. For each range of core numbers, the total mass surface density, Ed, of the cores are indicated by 
color-coded thin solid, dotted, dashed, long-dashed, dot-dashed, dot-long-dashed lines in order of increasing number (e.g. CI (solid) to C5 
(dot-dashed); C6 (solid) to C9 (long-dashed)). Squares mark the location of the best-fit cores. Heavier lines extend inward from near the 
squares show Ec,pl + ^cl,env5 i-e. based on the fitted power law density profile. 
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Fig. 16. — Same as Fig.[T5l now combining all 42 cores together (colored lines). The black symbols and lines show the masses of 31 actively 
star-forming core/clumps from Mueller et al. (2002), with triangles indicating the masses above a density threshold of nn ~ 3 x lO^cm"^ 
(we have scaled the masses by a factor 1.56 to be consistent with our adopted gas-to-dust mass ratio) and the squares indicating the masses 
inside the deconvolved source size. Note, the properties of the clouds on this inner scale are not directly resolved, but are inferred based 
on simple ID radiative transfer modeling. The IRDC cores/clumps overlap only with the lower-E range of the star-forming core/clump 
sample, perhaps indicating there is an evolutionary growth in core/clump density as star formation proceeds. 



TABLE 1 

Infrared Dark Cloud Sample^ 



Cloud Name 




b 


d i?eff 


e 


P.A. /,,o,obs^ /.,fore /J^^e ^(^^t) SgMF^ 


VBT09 
^SMF 


^SMF 


MBT09 








kpc pc 




° MJy/sr g/cm^ g/cm^ 


g/cm^ 


Mo 


Mo 



A (G018.82-00.28) 


18.822 


-0.285 


4.8 


10.4 


0.961 


74 


B (G019.27+00.07) 


19.271 


0.074 


2.4 


2.71 


0.977 


88 


C (G028.37+00.07) 


28.373 


0.076 


5.0 


15.4 


0.632 


78 


D (G028.53-00.25) 


28.531 


-0.251 


5.7 


16.9 


0.968 


60 


E (G028.67+00.13) 


28.677 


0.132 


5.1 


11.5 


0.960 


103 


F (G034.43+00.24) 


34.437 


0.245 


3.7^ 


3.50 


0.926 


79 


G (G034.77-00.55) 


34.771 


-0.557 


2.9 


3.06 


0.953 


95 


H (G035.39-00.33) 


35.395 


-0.336 


2.9 


9.69 


0.951 


59 


I (G038.95-00.47) 


38.952 


-0.475 


2.7 


3.73 


0.917 


64 


J (G053.11+00.05) 


53.116 


0.054 


1.8 


0.755 


0.583 


50 



93.8 


0.472 


0.209 


0.496 


0.0489 


0.0355 


17,700 


7,600 


85.1 


0.452 


0.075 


0.488 


0.0814 


0.0387 


2,200 


830 


89.9 


0.339 


0.266 


0.520 


0.0610 


0.0527 


45,000 


42,000 


71.5 


0.559 


0.327 


0.436 


0.0525 


0.0418 


53,400 


27,000 


96.7 


0.455 


0.276 


0.504 


0.0593 


0.0543 


25,200 


19,400 


48.4 


0.601 


0.193 


0.370 


0.0994 


0.0371 


4,460 


1,670 


43.2 


0.624 


0.140 


0.347 


0.0648 


0.0420 


2,010 


1,140 


45.6 


0.405 


0.142 


0.416 


0.0479 


0.0262 


13,340 


6,800 


42.1 


0.418 


0.141 


0.402 


0.0707 


0.0616 


2050 


1,490 


28.6 


0.509 


0.121 


0.328 


0.125 


0.0699 


259 


80 



^ Coordinate Names, Galactic coordinates, kinematic distances, effective radii (of equal area circles), eccentricities and position angles of fitted ellipses are from Simon et al. (2006). 
We then compare quantities derived in this paper with those from BT09. 

^ Mean intensity of the SMF background model (BT09) inside the Simon et al. (2006) ellipse. 

^ Areal average of those pixels for which values of SgMF > ^ derived. Estimates of a mean mass surface density based on Mgj^p and -Rgff are typically much smaller because of 
the regions inside the clouds ellipse with derived SgMF < (see BT09). 

The distance to IRDC F has recently been estimated from parallax of radio sources to be 1.56 ± 0.12 kpc (Kurayama et al. 2011), a factor of 0.42 smaller than the kinematic 
distance of Simon et al. (2006). For consistency with BT09, we retain the kinematic distance, but the mass of this cloud and its cores are likely to need to be reduced by a factor of 



TABLE 2 
Core Properties 



Core^ 


1 


b 


■Rcl = 








fep c 


"■H s cl 


"H s c 


Mc 




PL 




^cl 








"^H s cl 


"■H s c 




Mc 


■^c,PL 


^c PL 








pc 


g/cm2 


g/cm2 






lO^cm'-S 


lO^cm'-S 


Mq 


Mq 


g/cm2 


pc 


g/cm2 


g/cm2 






lO^cm'-S 


10^ ^ — 3 
cm 


JMq 


Mq 


Mq 


g/cm2 




18.78675 


-0.28592 


0.0962 


0.431 


0.316 




2.04 


2.47 


1.21 


37.9 


49.2 


0.353 


0.251 


0.319 


0.204 


1.38 


1.88 


0.713 


0.332 


303 


194 


204 


0.214 


A2^ 


18.80117 


-0.29625 


0.106 


0.356 


0.205 


1.30 


1.79 




0.838 


34.6 


35.9 


0.214 


0.251 


0.298 


0.149 


1.34 


1.42 


0.682 


0.353 


283 






0.162 


A3^ 


18.80750 


-0.30550 


0.109 


0.336 


0.249 


1.46 


2.06 


1.64 


0.767 


44.3 


46.0 


0.258 


0.698 


0.180 


0.0936 


1.48 


1.72 


0.135 


0.0594 


1320 


686 


687 


0.0937 


Bl*'"- 


19.28758 


0.08083 


0.114 


0.307 


0.216 


1.26 




1.61 


0.872 


42.1 


42.6 


0.218 


0.251 


0.240 


0.149 


1.28 


1.50 


0.570 


0.307 


228 


142 




0.148 


B2 


19.30758 


0.06625 


0.134 


0.222 


0.151 


1.10 


1.43 


1.08 


0.579 


40.7 


38.6 


0.144 


0.265 


0.185 


0.113 


1.12 


1.34 


0.452 


0.237 


195 


120 




0.109 


CI 


28.32450 


0.06655 


0.102 


0.383 


0.119 


0.900 


1.31 


2.72 


0.677 


18.6 


18.5 


0.118 


0.116 


0.382 


0.114 


0.900 


1.38 


2.40 


0.558 


77.8 


23.3 


23.4 


0.115 




28.34383 


0.06017 


0.0884 


0.510 


0.164 


1.20 


1.96 


3.59 


0.453 


19.2 


13.1 


0.111 


0.0872 


0.512 


0.165 


1.26 


1.42 


3.53 


1.71 


58.6 


18.9 


31.3 


0.274 


C3 


28.35217 


0.09450 


0.155 


0.165 


0.0686 


1.08 




0.711 


0.228 


24.9 


23.9 


0.0659 


0.145 


0.167 


0.0707 


1.08 




0.773 


0.244 


53.4 


22.5 


20.8 


0.0654 




28.35417 


0.07067 


0.0958 


0.434 


0.171 


1.09 


1.80 


2.99 


1.67 


49.5 


53.3 


0.386 


0.116 


0.421 


0.158 


1.22 


1.46 


2.24 


0.789 


85.8 


32.2 


34.8 


0.171 


C5 


28.35617 


0.05650 


0.104 


0.369 


0.0372 


0.781 


1.30 


2.72 


0.188 


6.04 


5.43 


0.0335 


0.0581 


0.403 


0.0512 


0.800 


1.70 


5.28 


0.390 


20.5 


2.60 


2.56 


0.0504 


C6 


28.36267 


0.05150 


0.109 


0.334 


0.177 


0.994 


1.82 


2.13 


0.718 


31.7 


34.3 


0.192 


0.581 


0.258 


0.103 


0.860 


1.26 


0.330 


0.102 


1317 


524 


501 


0.0985 


C7 


28.36433 


0.11950 


0.129 


0.238 


0.0556 


1.15 


1.45 


1.18 


0.718 


49.2 


43.3 


0.172 


0.0581 


0.277 


0.0727 


1.04 


1.58 


3.23 


0.629 


14.1 


3.70 


3.78 


0.0744 


C8 


28.38783 


0.03817 


0.111 


0.322 


0.202 


1.16 


1.58 


1.85 


0.895 


37.5 


37.5 


0.202 


0.290 


0.257 


0.137 


1.16 


1.46 


0.567 


0.245 


327 


175 


169 


0.133 


C9 


28.39950 


0.08217 


0.111 


0.321 


0.160 


1.02 


1.62 


1.99 


0.705 


29.8 


30.4 


0.163 


0.261 


0.277 


0.118 


1.02 




0.732 


0.267 


286 


121 


126 


0.122 


Dl 


28.52717 


-0.25033 


0.109 


0.337 


0.167 


0.956 


1.39 


2.17 


0.885 


29.7 


31.0 


0.178 


0.133 


0.325 


0.154 


0.960 


1.68 


1.71 


0.514 


86.0 


40.9 


39.9 


0.151 


D2 


28.53750 


-0.27650 


0.105 


0.361 


0.172 


0.873 


1.77 


2.53 


0.720 


28.6 


29.5 


0.0622 


0.332 


0.280 


0.0949 


0.860 


1.22 


0.623 


0.171 


464 


156 


153 


0.0930 


D3 


28.54150 


-0.23517 


0.110 


0.331 


0.0608 


0.788 


1.55 


2.29 


0.280 


11.0 


11.2 


0.243 


0.0663 


0.370 


0.0863 


0.720 


2.02 


4.38 


0.422 


24.4 


5.71 


5.46 


0.0825 


D4P 


28.53950 


-0.26950 


0.112 


0.318 


0.231 


0.772 


1.29 


2.18 


1.27 


43.5 


45.7 


0.262 


0.0663 


0.346 


0.254 


0.720 


1.32 


4.10 


2.10 


22.8 


16.8 


15.9 


0.241 




28.56533 


-0.22783 


0.0948 


0.443 


0.259 


0.880 


1.92 


3.43 


1.03 


35.1 


35.4 


0.0461 


0.531 


0.324 


0.145 


0.880 




0.445 


0.148 


1370 


614 


619 


0.146 


D6 


28.55550 


-0.23917 


0.126 


0.250 


0.0495 


0.586 


1.24 


1.66 


0.222 


11.9 


11.0 


0.188 


0.133 


0.249 


0.0486 


0.700 


1.62 


1.48 


0.168 


66.0 


12.8 


12.5 


0.0475 


D7^ 


28.56667 


-0.23300 


0.0993 


0.404 


0.168 


0.559 


1.28 


3.43 


1.12 


24.9 


27.9 


0.171 


0.265 


0.364 


0.133 


0.640 


1.54 


1.12 


0.247 


385 


141 


137 


0.130 


D8« 


28.57283 


-0.23267 


0.104 


0.369 


0.169 


0.900 


1.68 


2.57 


0.747 


27.6 


27.7 


0.0862 


0.497 


0.271 


0.0756 


0.840 


1.30 


0.408 


0.0856 


1010 


281 


269 


0.0724 


D9 


28.58667 


-0.22767 


0.126 


0.250 


0.0832 


0.834 


1.50 


1.49 


0.355 


19.9 


20.6 


0.0701 


0.232 


0.231 


0.0664 


0.840 


1.40 


0.743 


0.165 


187 


53.8 


56.1 


0.0693 


El 


28.64350 


0.13817 


0.115 


0.301 


0.0736 


0.947 


1.26 


1.86 


0.383 


14.6 


14.6 


0.0735 


0.0890 


0.313 


0.0826 


0.940 


1.36 


2.50 


0.561 


37.3 


9.85 


10.5 


0.0886 


E2 


28.64850 


0.12483 


0.127 


0.246 


0.164 


1.33 


1.66 


1.12 


0.678 


42.2 


45.1 


0.185 


0.267 


0.190 


0.109 


1.38 


1.72 


0.399 


0.177 


204 


117 


115 


0.107 




28.65883 


0.14350 


0.106 


0.352 


0.0983 


0.980 


1.78 


2.33 


1.35 




29.9 


0.176 


0.0890 


0.335 


0.0818 


1.48 


2.82 




0.175 


89.1 


21.8 


29.9 




Fl 


34.41950 


0.24583 


0.141 


0.199 


0.114 


1.36 


1.91 


0.806 


0.317 


34.4 


35.5 


0.118 


0.559 


0.139 


0.0549 


1.34 


1.60 


0.143 


0.0466 


658 


258 


253 


0.0537 


F2^P 


34.43517 


0.24217 


0.118 


0.284 


0.251 


1.36 


1.90 


1.37 


0.828 


53.0 


53.8 


0.255 


0.473 


0.184 


0.151 


1.32 


1.76 


0.227 


0.126 


622 


512 


468 


0.139 


P3sm 


34.44383 


0.24967 


0.127 


0.248 


0.178 


1.50 


2.01 


1.01 


0.513 


43.2 


46.2 




0.731 


0.140 


0.0718 


1.46 


1.76 


0.102 


0.0419 


1133 


578 


575 


0.0714 


F4 


34.45800 


0.25650 


0.128 


0.241 


0.122 


1.31 


1.98 




0.322 


30.4 


28.8 


0.116 


0.344 


0.184 


0.0664 


1.36 


1.66 


0.304 


0.0832 


329 




110 


0.0620 



34.73417 
34.78117 
34.78483 
35.47800 
35.48283 
35.48733 
35.48867 
35.49450 
35.52250 



0.56683 
0.56817 
0. 55750 
!1033 
D. 28700 
D. 29367 



0.377 
0.212 
0.232 



0.0723 
0.0674 
0.0910 



0.0785 
0.0605 
0.0166 



0.303 
0.289 
0.756 
0.849 
0.0872 
0.437 



0.0773 
0.0546 
0.0107 



0.392 
0.290 
0.245 



0.101 
0.0717 
0.123 



0.0943 
0.0730 
0.0337 



0.694 
0.366 
0.878 



0.725 
0.398 
0.258 
0.471 



42.0 
4.97 
37.8 
29.9 
36.2 



0.0945 
0.0718 
0.0323 



^ Core designation from BT09. ^ indicates core is saturated; P indicates core's central position has been moved out of the original core boundary from Rathborne et al. (2006); indicates the intensity minimum within the IRDC boundary occurs in this 
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TABLE 3 

Properties of IRDC Core and Clump Sample 



Mass ocale 


Property 


Mean 


Median 


Dispersion 


Mci(< r) = 60 M©^ 


= (pc) 


0.121 


0.114 


0.0238 




Eel (g cm-2) 


0.296 


0.318 


0.0952 




kp,c\ 


1.09 


1.10 


0.236 




nHscl (lO^cm-S) 


1.76 


1.85 


0.852 




' Mc(M0) 


30.4 


30.4 


12.7 




Ec (gcm-2) 


0.139 


0.160 


0.0738 






1.64 


1.67 


0.271 




s c (lO^cm-^) 


0.639 


0.750 


0.394 




_Mc,PL (Mo) 


31.1 


31.0 


13.5 




Sc,PL (gem 2) 


0.154 


0.171 


0.0899 


Best-fit PL+env.b 


Mel (Mo) 


279 


86.0 


392 




i?c (pc) 


0.222 


0.134 


0.187 




Eel (gcm-2) 


0.274 


0.277 


0.0925 






1.10 


1.12 


0.246 




nu,s,ci (105cm-3) 


1.76 


1.12 


1.63 




Me (Mo) 


128 


32.2 


185 




Ec (|cm-2) 


0.113 


0.103 


0.0535 






1.58 


1.56 


0.277 




nu s c (lO^cm-^) 


0.515 


0.354 


0.616 




Mc,PL (Mo) 


128 


34.8 


181 




Ee,PL (g Cm-2) 


0.121 


0.107 


0.0708 



^ Mass scale set by the enclosed mass in the E map being equal to 60 Mq . 
^ Mass scale set by the best-fit power law plus clump envelope model. 



